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The mass of the lowest spin-zero, strangeness- (— 2) flavor singlet state in the dibaryon sector has been calculated 
in quenched QCD on 16 3 x 32 and 24 3 x 32 lattices at f3 = 5.85 to study whether the energy of the proposed 
H dibaryon is near or below the AA threshold. Preliminary results indicate that finite lattice volume artifacts 
overestimate the binding, and that on the largest lattice mu is of the order of 100 MeV above the AA threshold. 



1. Introduction 

Based on physical bag model arguments and 
the magnetic hyperfine interaction, Jaffc sug- 
gested that the lowest bound state in the dibaryon 
sector would be a spin strangeness —2 SU(3) 
flavor singlet, and that the energy of this H 
dibaryon could be near or below the AA thresh- 
old [Q. Despite two decades of effort, experi- 
ments have neither found it nor excluded the pos- 
sibility of exotic multiquark systems (Q n Q ,n , for 
n + m > 3). 

If a six quark system exists as a bound state, a 
single gluon exchange will contribute 
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ij 



(oj • ffj)(A, ■ Xj)Mi : 
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to the mass jjj yielding Jaffe's original bag model 
estimate of m H = 2150 MeV, 81 MeV below 
the 2231 MeV AA threshold. Subsequent bag 
model predictions [|| yield the mass range 1.03 — 
2.3 GeV. Although early lattice calculations dis- 
cussed below were not definitive i the success 
of contemporary spectroscopy in quenched QCD 
motivated us to attempt a new lattice calculation 
of the H dibaryon. 

2. Lattice Parameters 

We have studied two lattice sizes, 16 3 x 32 and 
24 3 x 32 at quenched (3 = 5.85, corresponding to 
a lattice spacing a « 0.13(3) fm. Wilson fermion 
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propagators were calculated at seven values of the 
hopping parameter for the bare quark mass in 
the range 30 — 300 MeV using a point source at 
t = 0. For the fermions the boundary conditions 
were periodic in spatial directions and hard wall 
(Neumann) in the temporal direction. This choice 
allowed us to obtain a signal at large euclidian 
times which is crucial for the calculation. A total 
of 56 configurations at 16 3 x 32 and 40 at 24 3 x 32 
were used. 

Once the quark propagators were calculated, 
the two-point functions for the H as well as for 
assorted mesons (tt, p, <fi, K and K*) and hadrons 
(N, A and S) were constructed with the three- 
momentum p — at the sink. The lattice spacing 
is defined by ri^r, the degenerate u and d quark 
masses are defined by m n /m^ and the strange 
quark mass is defined by m^/uiN- These defini- 
tion were chosen to absorb as many of the errors 
arising from the quenched approximation as pos- 
sible. 

3. H Two-point Function 

While the H propagator could be expressed 
in terms of two-hadron states, this form is 
not computationally feasible due to numer- 
ous duplicate terms and mutual cancellations. 
Instead, we constructed the polynomial aris- 
ing from contracting J = 0, S = —2, 
B = 2 sources and sinks of the form <f>H = 
Cie abc e a ' b ' c ' eape^^xiuuddss)"^^, , and ma- 
nipulated it analytically to reduce the computa- 
tional complexity by a factor of 200. The final ex- 
pression was calculated by two independent mcth- 
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ods and symmetries of the result were checked. 

As a result, the machine time was divided 
evenly between the Dirac inverter and two-point 
function calculations for 16 combinations of heavy 
and light quark masses. 

4. Measurements 

In the absence of finite boundary effects, one 
may identify the lowest mass by writing the zero 
momentum two point function 

D(t) = J d 3 xe lpx D{0;x,t)\ p = 

= 5> n e-<"l^>< (2) 
|n) 

where n runs over all states with given quantum 
numbers and extracting the exponent with small- 
est mass. In our case, there is an additional cor- 
rection arising from the hard wall, so we approxi- 
mate D(t) in the region of interest by three terms: 

D(t) ra a ie - mit + a 2 e~ m2t + a 3 e- m3t , (3) 

and fit data in the window t\ < t < t%. In this 
case, mi represents contributions from the ex- 
cited states, to 2 is the ground state mass and 
TO3 accounts for the wall effects. This fit has 
been chosen over a simpler one mass fit because 
it is less sensitive to the window size and posi- 
tion. In addition, the resulting values of mi are 
less contaminated by higher excitations and wall 
effects. Typical results of the fit to the A and 
H on the 24 3 x 32 lattice are shown on figures [j] 
and H respectively. As one can see, the wall ef- 
fects on baryons are rather small, while, as ex- 
pected from the much more rapid fall off for the 
dibaryon, there is a larger effect on the H . We 
note that in all cases, eq. ([|) yielded a good fit to 
the linear interior region that was stable against 
reasonable changes in the {ti,t-z) window. Errors 
in the mass mi were determined by least squares 
analysis. Physical masses were determined by lin- 
ear extrapolation in the light and strange quark 
masses. 

5. Baryon Spectrum 

Since m u and rrid are fixed by the nucleon mass, 
in the non-strange sector, m p is a test both for 
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Figure 1. Lambda Propagator 
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Figure 2. Dibaryon Propagator 



finite size effects and quenching errors. On both 
lattice sizes we obtained comparable results, and, 
as expected ||, m p is about 10% lighter than the 
experimental value. 

In the strange sector, having set the strange 
quark mass by m^, the lattice results [and cor- 
responding experimental values] are: = 
990(6) MeV [1020 Me V], m K = 600(5) MeV 
[(m K + + m K o)/2 = 496 MeV], m K , = 
840(20) MeV [894 MeV] and m H = 1360(34) MeV 
[(m H o + m H -)/2 = 1318 MeV]. The errors are 
comparable with other quenched calculations. 
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Table 1 

Masses of the H and A in 24 3 x 32 lattice. 



ma 







H 


A 


0.1520 


0.1568 
0.1580 
0.1592 
0.1601 


2.091(09) 
1.988(10) 
1.858(14) 
1.754(20) 


1.0143(06) 
0.9602(08) 
0.9053(17) 
0.8586(24) 


0.1539 


0.1568 
U.lOoU 
0.1592 
0.1601 


2.021(08) 

1.799(15) 
1.680(22) 


0.9790(07) 

0.8681(20) 
0.8209(23) 


0.1554 


0.1568 
0.1580 
0.1592 
0.1601 


1.968(07) 
1.878(11) 
1.753(17) 
1.642(26) 


0.9501(08) 
0.8950(12) 
0.8371(23) 
0.7876(23) 


0.1568 


0.1568 
0.1580 
0.1592 
0.1601 


1.916(07) 
1.823(12) 
1.705(22) 
1.605(33) 


0.9221(10) 
0.8647(15) 
0.8060(26) 
0.7559(25) 



6. H mass 

For the lattice sizes we are able to treat, mn 
has strong finite lattice size effects. 

On the 16 3 x 32 lattice, the extrapolation to the 
physical m u and m s gives mn = 1950(60) MeV, 
just above the deuteron threshold. However, the 
2 fm lattice appears too small to rule out consid- 
erable finite size effects. 

Results from the larger lattice, 24 3 x 32 at 
the same value of (5 yield an unbound H at 
2340(20) MeV. The dibaryon remain heavier than 
two A's at all combinations of the hopping param- 
eters, as shown in table [j]. 

The sign and magnitude of the finite volume 
dependence indicate that interactions with peri- 
odic images are attractive, and significant at a 
separation of 2 fm. Since the attraction decreases 
with increasing volume and the H is already un- 
bound by the order of 110 MeV on 24 3 x 32 lat- 
tice, we believe the infinite volume extrapolation 
of quenched QCD is even more unbound. 



7. Comparison with Previous Work 

There were two previous attempts to calcu- 
late ma on the lattice a decade ago. Macken- 
zie and Tacker |3j did not see a bound dibaryon 
on a rather small lattice 6 x 6 x 12 x 18. Three 
years later Iwasaki, Yoshie and Tsuboi [Q used 
a larger lattice (16 4 x 48) and obtained a very 
strongly bound H, between 1450(250) MeV and 
1710(140) MeV, depending on the extrapolation. 
Since they used an unconventional gauge action 
and quark masses considerably heavier than those 
used in our calculation, direct comparison of nu- 
merical results is difficult. However, their result 
is consistent with the fact that interactions with 
images tend to overbind. 

8. Conclusion and Future Plans 

To strengthen our preliminary conclusion that 
the H is unbound in quenched QCD, we intend 
to recalculate the small lattice propagators at the 
same hopping parameters at the large lattice, so 
that more direct study of the finite size effects 
can be made. We also will calculate the two- 
point functions for dibaryons with other quantum 
numbers. In addition, we will extrapolate hadron 
masses in and mx rather than in the quark 
masses. 
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